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Abstract 

This paper examines the properties of the homentropic Euler equa- 
tions when the characteristics of the equations have been spatially aver- 
aged. The new equations are referred to as the characteristically averaged 
homentropic Euler (CAHE) equations. An existence and uniqueness proof 
for the modified equations is given. The speed of shocks for the CAHE 
equations are determined. The Riemann problem is examined and a gen- 
eral form of the solutions is presented. Finally, numerically simulations 
on the homentropic Euler and CAHE equations are conducted and the 
behaviors of the two sets of equations are compared. 



1 Introduction 

This manuscript examines a modification of the homentropic Euler equations 
where the characteristics of the equations have been spatially averaged. The 
equations examined are 

- ^ 

Pt + upx + p-Ux = (la) 
_ ad 

ut^uux^ Px=0 (lb) 

P 

u = g ^ u (Ic) 
a = g ^ a (Id) 

with = 7/)^~^. These are derived in section [2j The inspiration for these equa- 
tions came from previous work done on a similarly modified Burgers equation. 

The Burgers equation, Ut + uux = 0, is considered a simplistic model of 
compressible flow. Multiple independent investigations have been made into a 
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modified Burgers equation with an averaged convective velocity [H [21 |3l [H |5l [6] , 



Ut + uux = (2a) 
u = g'^ (2b) 

= -9{-). (2c) 
a a 

where g is an averaging kernel. A brief summary of the primary results are as 
follows. 

1. Solutions to Equations ([2| exist and are unique [11 12]. 

2. When the initial conditions are C^, the solution remains for all time 

i- 

3. For bell shaped initial conditions it was proven that as the averaging 
approaches zero (a 0), the solutions to Equation ([2| converge to the 
entropy solution of inviscid Burgers Equation [5] . 

4. There are a set of discontinuous initial conditions where solutions to Equa- 
tion ([2| will converge to a non-entropic solution, but these solutions are 
unstable [5, 3 . 

5. It is conjectured that solutions to Equation ([2| will converge to the entropy 
solution of inviscid Burgers Equation for all continuous initial conditions, 
and thus the scheme 

Ut + uUx = (3) 
u = (4) 
ix(x,0) ^uq{x), (5) 

will converge to the entropy solution for any bounded initial condition uq 

The work on the regularization of the Burgers equation is inspired by and 
related to work done on the LANS-a equations [1 [1 [101 [TH [El [13] . These 
equations also employ an averaged velocity in the nonlinear term and have been 
successful in modeling some turbulent incompressible flows. 

It is thought that a similar regularization could be accomplished for the 
equations that describe compressible flow. Encouraged by the results for Burg- 
ers equation, the next step is to attempt to introduce averaging into the one- 
dimensional homentropic Euler equations, a simplified version of the full Euler 
equations, where pressure is purely a function of density. There have been 
several attempts at such a regularization. 
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Using a Lagrangian averaging technique Bhat and Fetecau [14] derived the 
fohowing equations 




(6) 



(7) 



(8) 
(9) 



While the solutions to the system remained smooth and contained much struc- 
ture it was found that the equations were " not well-suited for the approximation 
of shock solutions of the compressible Euler equations." 

Another attempt by Bhat, Fetecau, and Goodman used a Leray-type aver- 
aging [l5 leading to the equations 



with p = np^ . They then showed that weakly nonlinear geometrical optics 
(WNGO) asymptotic theory predicts the equations will have global smooth 
solutions for 7 = 1 and form shocks in finite time for 7 7^ 1. 

Additionally in 2005, H. S. Bhat et. al. [16 applied the Lagrangian aver- 
aging approach to the full compressible Euler equations. Their approach was 
successful in that a set of Lagrangian Averaged Euler (LAE-a) equations were 
derived. However, the equations seemed so intractable that they seemed im- 
practical for real world applications. 

Inspired by the existence uniqueness proofs from the averaged Burgers equa- 
tions found in 01 E] , we average the characteristics of the homentropic Euler 
equations to derive what we term as the characteristically averaged homentropic 
Euler (CAHE) equations ([T]). It is the properties of these equations that this 
paper examines. 

The following section follows the derivation and present the final equations. 
The existence and uniqueness of solutions to the equations are then proven in 
section [3j Sections [4] and [5] examine the speed of the shocks and solutions to 
the Riemann problem. Numerical simulations and their comparison to those of 
the homentropic Euler equations are discussed in sections [6] and jT] The results 
are then briefly summarized in the concluding remarks. 

2 Derivation of the equations 
2.1 The homentropic Euler equations 

We begin the process of deriving the CAHE equations by starting with the 
homentropic Euler equations. There are two equations to the homentropic Euler 



Pt + upx + pUx = 

Px 

Ut + uux H = 

P 



(10) 
(11) 



u = u — a^u. 



(12) 
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equations. Conservation of mass and conservation of momentum. Pressure is 
expressed purely density raised to the power of 7. 



Pt + {pu)x = 

{pu)t + {pUU^ p^):r: = 

The equations are then written in primitive variable form 



0, 



" p 




u 


P 




' P 




u 


t 


P 


U 




u 


X 



(13a) 
(13b) 



(14) 



with o? = 7p^~^. 

To get a clear view of the characteristics involved the equations are then 
diagonalized to obtain the equations 



where 







u-\- a 




+ 

t 







= u -\- 




V 


= u — 





u — a 



2a 

2a 
7^* 







(15) 

(16a) 
(16b) 



The variables are commonly know as Riemann invariants. From Equation 
( 15 ) it is easy to see that the quantity is convected at speed u -\- a and v~ 
speed u — a. Thus along the characteristics + a, will remain constant. 

Shocks will form when characteristics intersect. To prevent this from hap- 
pening the characteristics are spatially averaged. This averaging is conducted 
by convoluting the variable to be averaged with an averaging kernel, ^, and is 
represented by a bar above the variable. For example, the averaged velocity 
would be expressed 

u = g ^u. (17) 



2.2 The averaging kernel 

Several assumptions on the averaging kernel are made at this point. The kernel 
is assumed to be even, for isotropic purposes. For the purpose of Theorems 



3.1 and 3.2 the kernel and its first derivative are assumed to be integrable. Of 
special interest is the Helmholtz filter which is defined as 



and thus has an averaging kernel of 

g{x) = iexp(-|x|). 



(18) 



(19) 
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For all numerical simulations found in Sections [Sj |6j and [7| the Helmholtz filter 
is used for its convenient inversion techniques. 

Furthermore the kernel will be equipped with a parameter a which will 
control the amount of averaging. If g is the averaging kernel, then a is introduced 
as 

5" = -g (-) ■ (20) 

Thus as a ^ the averaging kernel approaches the Dirac delta function. 



2.3 Obtaining the CAHE equations 

Using the averaging kernel discussed above the characteristics of the homen- 
tropic Euler equations are averaged to obtain the CAHE equations in diagonal- 
ized form. 

r 1 r I 7^ n 1 r 1 

(21) 





+ 




t 



u -\- a 
u — a 



0. 



Change the equations back into the primitive variable form and we get the 
equations 



Pt + UPx + p-Ux = 

a 



Ut + UUx + 



aa 







P 

u = g ^ u 



a 



g^a 



with = jp^~^ which were seen earlier as Equations ([I]). 

These are now the equations that will be examined for the rest of this paper 
and are referred to as the characteristically averaged homentropic Euler (CAHE) 
equations. We begin by proving that Equations ([T]) have one and only one 
solution. 



3 Existence and uniqueness theorem 

A critically important property for the Equations ([T]) to have is that a solu- 
tion exists. This section addresses this problem by presenting a proof for the 
existence and uniqueness of solutions. This proof uses a method of character- 
istics approach and is similar to the existence uniqueness proof of Convectively 
Filtered Burgers equation found in [4 

An outline of the proof is as follows. When the equations are cast into their 
characteristic form, it is clear to see that the fluid velocity and the speed of 
sound remain bounded. Since those speeds remain bounded, the first derivative 
of the averaged speeds will remain bounded as well. The characteristics are 
governed by the averaged speeds and with the first derivatives bounded, the 
characteristics may grow closer, but will never intersect. Thus the solution can 
be fully realized by the characteristics and the initial conditions. 
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Theorem 3.1 Let g{x.) G ly-^'-^(R) and uo{x)^ao{x) G C'^{R), then there exists 
a unique global solution u{x , t) ^ a{x , t) G C-^(R,M) to the initial value problem 



Pt + upx + p-Ux = 
a 



aa 

Ut + UUx H px 

P 







u = g ^ u 
a = g ^ a 
u{x, 0) — Uq 
a(x, 0) = ao 



(23a) 

(23b) 

(23c) 
(23d) 
(23e) 
(23f) 



with = 7P 



.7-1 



Proof Express Equations ( |23a| ) and ( |23b| ) with matrices as 





" p ' 




u 

aa 


P- 1 




' P ' 


u 






U 




U 




t 


- P 









By diagonahzing the matrix these equations can be rewritten as 

n + 





+ 




t 



u-\- a 
u — a 



0. 



(24) 



(25) 



It is here that we shift our perspective to a method of characteristics type 
view. Associate the maps (j)^ with the characteristics of v^. Thus we have 



9 + 



and 



d_ 



d 



along a (f)'^ characteristic 
along a characteristic . 



From this we can obtain the estimates that 



\v^{x,t)\ 



(26) 
(27) 

(28) 
(29) 

(30) 



It the mappings <p^ have continuously differentiable inverses, ip^, then Equa- 
tion (25) has the solution 



(31) 
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Sufficient conditions for such inverses to uniquely exist is if the Jacobians of 
(j)^ are non-zero for ah positions and time. Thus if J{(l)^) 7^ 0, Equation (25) is 
uniquely solved by (31). 

For the next section of the proof we will be dealing with The results 
for follow in precisely the same manner. Since we are dealing with ID the 
Jacobian of is essentially It is clear from Equation (26) to see that the 
time derivative of ^+ is 



Thus we see that 



dt 

= ^J(O) exp 



Thus (bt will remain non-zero if 



Jq Ux + ax dt 



a^r dt 



< 00. 



(32) 
(33) 



First we will show that u^a G . We see that 



\2u\ 



2a 



< 



7-1 
.,+ 11 



-\- u -\- u ■ 



2a 



7-1 



(34) 
(35) 



Looking at Equation (30) and noting that ||'^^||/^oo and H'^^H^^oo are bounded 
for all time, then ||i^||^oo is bounded for all time. Similarly one can bound the 
quantity ||a||^oo for all time. 

Given that g{x) G W^'^(R), there exists M G M, such that 



d_ 



< M < 00. 



(36) 



Knowing that -^g G and u G L^, we know that -^u exists and that 

d d 



—-g * u. 
ox ox 



Similarly 



d d 
-^CL = —-g * a. 
ox ox 

Using Young's inequality we can bound the derivatives Ux and clx. We get 

(37) 
(38) 



d _ 




d 






dx 


< 


dx^ 




d _ 




d 






dx 


< 


dx^ 





This leads directly to the bound 



Jo 



ttx dt 



+ \\a\\L^)t 



(39) 
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Thus for finite time, tiie Jacobian of 0^ remains uniquely invertible, with a 



continuously differentiable inverse and thus (31) is a unique C^(M^) solution to 



(25). A unique C^(M^) solution to (23) follows accordingly. 



Theorem 3.2 Let g{x.) e W^'^(R) and uo{x),ao{x) e L^(R), then there exists 
a unique global solution u{x ^ t) ^ a{x ^ t) G L'^(R,R) to the initial value problem 



a 

pt + UPx + p-Ux = 

a 


(40a) 


ad 

ut + uux H px = 

P 


(40b) 


u = g ^ u 
d = g ^ a 
u{x, 0) = Uq 
a(x, 0) = ao 


(40c) 
(40d) 
(40e) 
(40f) 



Proof The proof is the same as for Theorem 3.1 (j) still have unique con- 



tinuously differentiable inverses, and the solution remains in the same form as 



Equation (31), but now lacks continuity due to the initial conditions. 



It should be noted that different averagings can be used for u and d. The 
existence and uniqueness proof will hold as long as a bound on the first deriva- 
tives of u and d remains constant. This is noted as the type of averaging used 
will affect shock speed as noted in Section |4] 

4 Shock speeds 

One of the consequences of the previous sections is that for continuous initial 
conditions, the solution will remain continuous. Also for initial conditions with 
a discontinuity, that discontinuity will remain. There are no colliding charac- 
teristics. Thus the "shock" is really a discontinuity that is being convected. 
If there is a discontinuity in Vq{x) then that discontinuity will travel at speed 
'u(x*)+a(x*), where is the location of the discontinuity. Similarly, if there is a 
discontinuity in Vq (x) then that discontinuity will travel at speed u{x^) — a(x*). 

Consider the case where there is a single jump discontinuity in but 
otherwise constant. 

v+{x) = ( ""tn (41a) 
v-{x) =C. (41b) 
This will lead to a jump discontinuity in u and in a. 



With u and a defined as in Equations (23) and the filter g being even and 



J g = 1 then the speed of the discontinuity in Equations (41) will be 



s = — ^ + — ^ — (42) 
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where ui^Ur^ai and are the hmiting values of u and a on the left and right 
side of the discontinuity respectively. 

Similarly, if the jump discontinuity existed only in v~ then the speed of the 
discontinuity would be 

ui^Ur ai^ar 
'=—2 (4=^) 

For multiple discontinuities or for nonconstant values around the discontinu- 
ities, the values of u and a cannot generally be analytically written . However, 
for small values of a the filter will not "see" as far and the values of u and 
a will be roughly the average of the values just to the left and right of the 
discontinuity. 



4.1 Alternative averagings and shock speeds 

In the above calculation u and a were defined using the definition 

(•) = <?*(•) (44) 

This is not, however, the only way to conduct averaging. One possible alterna- 
tive is using Favre averaging, a density weighted average, where 

(•) = (45) 

This alternative will not affect the existence uniqueness theorems established 
in section |3] It will, however, change the shock speeds established earlier. 
If Equation (45) is used, the speed of a discontinuity in v+ will travel at 

PlUi^PrUr , pm^Prttr 

S= \ . (46) 

Pi ^ Pr Pl^ Pr 

Another possibility is to apply a spatial average to the conserved quantities, 
density and momentum, and then compute the characteristic speeds from the 
the averaged conserved quantities, 

u = ^ (47) 
P 

~a'=j{pr-\ (48) 

If the speed of the characteristics were defined using u and a, then the speed of 
the shocks would be 

^^Pim+PrUr ^^(pi+PrV-^ ^ (49) 
Pi ^ Pr \ 2 J 

Thus how one chooses to define u and a will affect the speed of the discontinu- 
ities. 
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5 Riemann solutions 



This section builds upon the results of section [4] and establishes properties of the 
solutions of the CAHE equations for the Riemann problem. It then continues on 
to examine solutions of the CAHE equations if, instead of the Riemann problem, 
the initial conditions are slightly perturbed. 



5.1 Riemann problem 



Looking at Equation (|3l|, it is clear that any discontinuity that exists in the 
initial condition remains in the solutions for all time. Furthermore, no additional 
discontinuities will form. 

For the Riemann problem, at time t = a single discontinuity will exist at 
the origin: 



-{ 



a{x,0) = { "^JJJ (51) 



which can also be rewritten in terms of the variables and v . 



These discontinuities will travel at the speeds discussed in section |4] We 
denoted the locations of the discontinuities in v+ and v~ as x~ and The 
speeds of these discontinuities are determined by 

d 

—x^ = u{x^)±a{x^). (54) 

Clearly the discontinuity found in will travel at a faster speed than the dis- 
continuity in v~ . Thus the solution will consist of three different areas separated 
by the locations of x~ and x+. The solution is then 

v+{x,t) = \ ^1,fj^^ (55) 

-/ 4\ { v7 X < x~(t) 
V {x,t) = < L J ^ (56) 

^ ^ [ X [t) < X ^ ^ 

If expressed in primitive variables the solution is 

X < x~{t) 

u{x,0)={ Um x~{t) < X < x^{t) (57a) 

X~^{t) < X 

X < x~{t) 

a{x,0) = { am x~{t) < X < x^{t) (57b) 

X~^{t) < X 
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where Um = ^^"^ and = ^^"^ • 

Thus in general, the Riemann problem will produce the solutions presented 



in Equations (57) where two discontinuities will propagate dependent upon the 
speeds of the averaged velocity and speed of sound. If the initial conditions hap- 
pen to be chosen, such that or v~ are constant, then there will exist only one 
traveling discontinuity. The following subsection show numerical simulations of 
a typical Riemann solution. 

5.1.1 Riemann solution numerics 

In section [6| a numerical scheme is described that numerically simulates the 



behavior of the CAHE equations. However in section 5.2 it is demonstrated that 
some of the solutions to the Riemann problem are unstable and thus cannot be 
capture by the numerical scheme. Thus this numerical scheme was developed. 
Since the Riemann problem solutions are of the form presented in Equations 



(57), finding solutions, even unstable ones, to the Riemann problem reduces to 
tracking the locations of the discontinuities in time. The following numerical 
scheme does precisely this. 



To find the solution (57), the values ui^rn,r and ai^rn,r are all known, one 
simply has to determine the location of x^. This is done by using Equation 



(54). A typical iteration is as follows. 

1. The values of u{x^ti) and a{x^ti) are known on a uniform grid. 

2. The Helmholtz operator is then numerically inverted to give the values of 
u and a on the uniform grid. These values are then interpolated to give 
the values of ± a at . 

3. The positions of the x^, are then advanced in time by the values of u±a 
on x^ calculated in the previous step. 



4. The values of u{x, ^2) and a(x, ^2) are now known using Equation (57) and 
current positions of x^. 

The inversion of the Helmholtz operator is described in section [6] 

To examine a typical case, we will simulate the Riemann solution with the 
initial conditions 



uix,0) = [l (58a) 



«(^,0) = | I (58b) 



This set of initial conditions will be referred to as Example la. This is a standard 
shock tube problem where the velocity is zero initially and there is a jump in 
pressure. The initial conditions can be seen in Figure [l] With these values 
there are discontinuities in variables and v~ . These discontinuities are seen to 
propagate to the right and left respectively as seen in Figure [2] These numerical 
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simulations were conducted with 2^^ grid points on the domain [—1,1] with 
a = 0.02. 



2 

;s 1 - C3 



-1 





X 

(a) 





X 
(b) 



Figure 1: The initial conditions for the Riemann problem, (a) The initial ve- 
locity profile is 7x = 0. (b) The initial speed of sound profile is a simple jump 
discontinuity, indicating a higher pressure on the right. 



2- 
;s 1 - 






X 

(a) 



C3 





X 

(b) 



Figure 2: The solution to the Riemann problem with initial conditions Equation 



(58). There are two shocks found in both the velocity and speed of sound. This 
is the solution at time t = 0.2. 



5.2 The perturbed Riemann problem 

While the solution to the Riemann problem generally leads to one or two trav- 
eling discontinuities, we have found that some of these traveling discontinuities 
are unstable. 

There have been examinations of the Burgers equation with the character- 
istics averaged [U [H |5j [2j |3] . In two of these papers it was seen that certain 
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traveling discontinuities in these equations were unstable and if perturbed would 
become expansion waves [5l [3] . We have found similar behavior for the C AHE 
equations. 

Consider Example la as seen in Figure [2] The position of the leftmost shock 
is denoted x~ . The value of — a is greater to the right oi x~ than it is to 
the left. So clearly — a is increasing across the discontinuity at x~ . Thus 
the v~ characteristics are diverging there. This is typically indicative of an 
expansion wave. However, the spreading of this discontinuity remains simply a 
discontinuity. 

Suppose instead that instead of a strict discontinuity, the function was a 
continuous function, albeit extremely steep. The diverging characteristics would 
then spread this extremely steep gradient and gradually make it less steep, 
forming an expansion wave. Thus, this traveling discontinuity is unstable, and 
any amount of smoothing would lead to an expansion wave as opposed to a 
shock. 

Using numerics discussed in section |6| we made runs using the same ini- 
tial conditions as in Example la except the initial conditions were initially 
"smoothed", using the same Helmholtz filter as employed in the averaging of 
the characteristics. We will refer to this as Example lb. This smoothing of 
the initial conditions is similar to the approach used in [5] . Figure [s] shows 
the initial conditions where the smoothing is readily apparent in the a variable. 
Figure |4] shows the solution at time t = 0.2. The right shock remains unchanged 
from the unperturbed solution seen in Figure [2] However, an expansion wave 
is clearly seen on the left in the perturbed solution where previously there was 
a discontinuity. These simulations were run with a resolution of 2^^ grid points 
and a = 0.02. 



2- 
;s 1 - 




C3 



-1 





X 

(a) 





X 
(b) 



Figure 3: The initial conditions for the perturbed Riemann problem, (a) The 
initial velocity profile is u = 0. (b) The initial speed of sound profile is a 
smoothed pressure jump. 



To summarize, the solutions to the Riemann problem take the form of Equa- 



tions (57) and produces one or two traveling discontinuities, some of which may 
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3 



3 






X 

(a) 





X 

(b) 



Figure 4: The solution to the perturbed Riemann problem. The initial condi- 
tions are Equation (58) that have been filtered with the Helmholtz filter. There 
is clearly a shock and an expansion wave. This is the solution at time t = 0.2. 



be unstable. If instead the discontinuities in the initial conditions are replace 
with steep, but continuous gradients, then the results will begin to include ex- 
pansion waves in addition to the shocks. 



6 Numerics 



In section 5.1.1 a numerical method is presented that will solve the Riemann 
problem. Here a numerical method is established that can be applied to more 
general initial conditions. In section [7| we compare the behavior of the CAHE 
equations to the homentropic Euler equations. To accomplish this we use a nu- 
merical technique based on the characteristic structure of the CAHE equations. 
This technique is similar in nature to that used by Bhat and Fetecau [3 . Since 
the equations we are investigating are inspired by the use of averaged charac- 
teristics, our numerical method will track the characteristics as they evolve in 
time. 

Our numerical method is interested in solving Equations ([T]) with initial 
conditions 

v+(x,0) =v^{x) (59) 
v-{x,0) = Vq{x). (60) 

Since we cannot operate on an infinite domain we restrict the problem to the 
domain [a, b] and assume that on the boundaries 

^v+ix,t)=0 (61) 
^v-{x,t) = 0. (62) 



14 



From Equations (26) and (28) we know that the value of does not change 
along its characteristics, thus 



v^{^^{X))=vt{X), (63) 

where X is the label of the characteristic. 

Let X be the equally spaced grid points on the interval [—a, a]. Let be the 
positions of the characteristics Xi. These vector can also be considered as the 
values of ^ or as a dynamic grid. At time t = 0, x = = X~. Notationally 
let (x) be a vector containing the values of v"^ evaluated at points x and let 
v^(X^) be evaluated at the locations of the characteristics X^. 

The positions of the characteristics are iterated in time, thus giving us the 
solution. A typical iteration proceeds as follows. 

L The values ofv^{X^) are known and remain constant through time as a 
result of Equation ( [63| ). Using cubic splines these values are interpolated 

to 'U^(x). 

2. From 'u^(x), it is a simple task of finding the values of u and a on the 
grid X. The Helmholtz operator is then numerically inverted to give the 
values of u and a on x. These values are then interpolated back to give 
the values of ± a on X^. 

3. The positions of the characteristics, X^, are then iterated in time by the 
values of ± a on X^ calculated in the previous step. 

In step 2, the Helmholtz operator is numerically inverted. To give a brief 
explanation of how the operator is inverted first consider a finite difference 
approximation of a second derivative on a equispaced grid. 

(64) 

Thus the Helmholtz operator (1 — <^^^) can be represented as a quintdiagonal 
matrix. This is easily inverted allowing the calculation of ± a from ± a on 
an equispaced grid. 

A benefit to this method is that plotting the vectors X^ versus time gives 
a graph of the characteristic plane. 

7 Comparisons with homentropic Euler equa- 
tions 

The CAHE equations were based on the homentropic Euler equations so natu- 
rally we compare the behavior of the two. In particular this section examines 
two examples chosen specifically to compare and contrast the equations behav- 
ior. Example 2 was chosen so that the CAHE equations would have a single 
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traveling shock. The initial conditions for Example 2 are 



«(^,0) = | I (65b) 



Example 3 was chosen so that the homentropic Euler equations would have a 
single traveling shock. The initial conditions for Example 3 are 

Hx,0) = [ %.37440483... x > H 

«(^,0) = | I (66b) 



Examples 2a and 3a are the simulation of the homentropic Euler equations and 
2b and 3b address the CAHE equations. In both examples the two sets of 
equations are found to behave significantly differently. This is attributed to the 
the fact that the CAHE equations will not produce new characteristics as shown 
by the proof of Theorem |3.1[ while the homentropic Euler equations do produce 
new characteristics. 



7.1 Numerics for the homentropic Euler equations 

Numerical simulations of the homentropic Euler equations require a separate 
method than that described in section |6] because characteristics in the homen- 
tropic Euler equations collide and are created. For the numerical simulations 
of the homentropic Euler equations the Richtmyer method, a well-established 
if low order method, was utilized as described by [T^. This method is second 
order finite difference scheme and employs an artificial viscosity. Clearly this 
is not the most optimal numerical scheme for the Euler equations, but run at 
a sufficiently high resolution it will suit our purposes. For reference there were 
2i2 g^^^^ points on a [—1,1] domain. This method requires an artificial viscosity 
for stability when examining the Riemann problem. Several different values of u 
were tested to see that the value did not significantly affect the solutions on the 
time interval examined. For the numerical simulations shown here, the artificial 
viscosity was set at = 0.08. 



7.2 Example 2 

Example 2 ( [65| ) was chosen so that the CAHE equations would have a single 
traveling shock. Notice that with these initial conditions there is a discontinuity 
in the variable but v~ is constant. Since v~ is constant in the beginning it 
should remain constant for all time. This is a consequence of the fact that no 
characteristics are created or destroyed as time progresses. 
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7.2.1 Example 2b, the CAHE equations 



Figure |5] shows the simulation for Example 2b which examine the CAHE equa- 
tions. The simulation was conducted with a resolution of 2^^ and a = 0.02. 
Figures andjsb show the single discontinuity progressing to the right as ex- 
pected. Figures [sj) and[5]i show the values of and v~ . The discontinuity in 
v+ has traveled to the right and as expected v~ remains constant. 

Figures and [sf show the paths of the and v~ characteristics respec- 
tively. Figure shows that the characteristics of are moving towards each 
other and as they near the shock are bent toward each other. They do not, 
however, intersect. Figure |5]F shows that the characteristics of v~ pass through 
the shock and upon doing so change their speed to match the characteristics on 
the other side. 



7.2.2 Example 2a, the homentropic Euler equations 

Now we examine the behavior of homentropic Euler equations for the same 
initial conditions. Figure |6] shows the simulation for Example 2a. The simulation 
was conducted with a resolution of 2^^. 

Figures and [BfD show a shock progressing to the right as expected. In 
addition there is an expansion wave also occurring. Figures [gJ: and[6]i show the 
values of v+ and v~ . Of primary notice is that v~ is no longer a constant but 
has attained new value. This is a significant departure from the behavior of the 
CAHE equations. Figures and |6f are not graphs of the actual simulation. 
The Richtmyer method does not lend itself to characteristic graphs as does the 
method used for the CAHE equations. Instead they are sketches that depict 
the behavior of the giving simulation. 

The characteristics in change speed as they travel through the expansion 
wave. The characteristics intersect causing the shock. Thus at the shock the 
characteristics are destroyed. 

For the v~ characteristics, at time t = 0, the characteristics to the left of 
the discontinuity have a speed greater than those to the right, but both are less 
than the speed of the shock. The result is that the characteristics on the right 
will be absorb into the shock and between the characteristics on the left and the 
shock there will be a gap devoid of characteristics. This gap is filled with an 
expansion wave created at t = and new characteristics that are created at the 
location of the shock continually as time progresses. The values of these new 
characteristics can be determined with the Rankine-Hugoniot jump conditions. 
This creation of characteristics again is a behavior that the CAHE equations 
does not demonstrate. 



7.3 Example 3 

The previous example was chosen so that the CAHE equations would have a sin- 
gle shock. In contrast. Example 3 (66) was chosen so that the homentropic Euler 
equations would have a single traveling shock. Again a significant departure in 
behavior will be noticed. 
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(a) (b) 



18r -2r 

16- 




Figure 5: Example 2b at time t = 0.1. This example was chosen so that the 
CAHE equations would have a single traveling shock which is clearly visible in 
the velocity and speed of sound in (a) and (b). (c) and (d) show graphs of the 
Riemann invariants and v~ . Note that v~ is constant as it was in the initial 
conditions, (e) and (f) show the and v~ characteristics respectively. Notice 
that the characteristics are converging to the shock but never intersect, while 
the v~ characteristics pass through the shock and change speeds as they do so. 



7.3.1 Example 3a, the homentropic Euler equations 

First we examine the behavior of homentropic Euler equations for Example 3a. 
Figure [7| shows the simulation for Example 3a. The simulation was conducted 
with a resolution of 2^^. 

Figures [7^ and[7]3 show a single shock. The shock is close to stationary but 
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(a) (b) 
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(e) (f) 

Figure 6: Example 2a at time t = 0.1. While the CAHE equations form a single 
shock the homentropic Euler equations clearly form a shock and expansion wave, 
(a) and (b) show the shocks and expansion wave clearly in the velocity and 
speed of sound, (c) and (d) show the Riemann invariants and v~ . Notice 
that while v~ began as a constant it is no longer. This is due to the creation of 
new characteristics at the shock, (e) and (f) show sketches of the characteristics 
for this example, (e) shows that the characteristics are being absorbed by 
the shock, (f) shows the v~ characteristics being created at the shock. The 
dotted line represents the expansion wave. 



is progressing to the right. Figures [t]^ and [7]i show the values of and v~ . 
The discontinuity can be seen in both Riemann invariants. This is noticeably 
different than the behavior of the CAHE equations where a single traveling 
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shock will appear in only one of the invariants. 

As in the previous examples Figures [7^ and [Tf are not graphs of the actual 
simulation, but sketches that depict the behavior of the giving simulation. The 

characteristics to the left of the discontinuity travel faster than those to the 
right with the speed of the shock found in between those speeds. Thus the char- 
acteristics collide and cause the shock. The v~ characteristics to the left of the 
discontinuity travel faster than those to the right, however both are slower than 
the speed of the shock. Thus there is an area devoid of characteristics between 
the left characteristics and the shock. This is filled with new characteristics that 
originate from the shock. 

7.3.2 Example 3b, the CAHE equations 

Now we examine the CAHE equations for Example 3b. Notice that with these 
initial conditions there is a discontinuity in the both variables and v~ . Thus 
it is to be expected that there will be two phenomenon, either expansion waves 
or shocks. In this case, both will be shocks. 

Figure [5] shows the simulation for Example 3b. The simulation was con- 
ducted with a resolution of 2^^ and a = 0.02. 

Figures [5^ and|5]3 clearly show two distinct shocks, close together, progress- 
ing to the left. These shocks can also be seen in Figures (Sj) and |5]i in the 
variables and v~ . 

Figures and [sJF show the paths of the and v~ characteristics respec- 
tively. The line across the center of the graph represents the remapping of char- 
acteristics as they had grown too close together for convenient computation. 
For both sets of characteristics there is a convergence of the characteristics to- 
wards the shock. This demonstrates that it is a stable discontinuity and that a 
perturbation will not turn these into expansion waves. Again the characteristics 
are bent towards each other, but never intersect. 

7.4 Comparisons between the homentropic Euler and the 
CAHE equations 

From these examples and previous sections we can begin to draw comparisons 
between the homentropic Euler equations and the CAHE equations. 

Both sets of equations are capable of forming shocks and expansion waves 
from the Riemann problem. For homentropic Euler equations it is well estab- 
lished that for the Riemann problem there will often be both a shock and an 
expansion wave formed. With the CAHE equations, introduced in this paper, 
there will always be two shocks formed. However, if the initial conditions are 
slightly smoother, some of those discontinuities will be found to be unstable and 
result in an expansion wave. Thus both behaviors can be said to be found in 
both equations. 

In terms of speed of the shocks, we have found there to be differences. In 
the homentropic Euler equations, the speed of the shocks are determined by the 
Rankine-Hugoniot conditions, in order to preserve mass and momentum. With 
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(a) 



(b) 




(e) (f) 

Figure 7: Example 3a at time t = 0.2. This example was chosen to have a single 
shock for the homentropic Euler equations which is clearly visible in (a) and 
(b). The shock is moving to the right but clearly not as quickly as previous 
examples, (c) and (d) show the Riemann invariants and v~ . (e) and (f) are 
sketches of the characteristics for and v~ . Again in (f) v~ characteristics 
are being created at the shock. 



the CAHE equations the speed of a discontinuity is determined by the speed 
of the averaged characteristics at the location of the discontinuity. In section 
|4j it can be seen that with the averaging chosen for this paper that the shock 
speed differs from those of the homentropic Euler equations. From this it is also 
clear that while the homentropic Euler equations may be formally regained from 
the CAHE equations by letting a ^ 0, the solution will not converge to weak 
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(e) (f) 

Figure 8: Example 3b at time t = 0.2 (a) and (b) show that there are two 
distinct left traveling shocks in the velocity and speed of sound, (c) and (d) 
show that these shocks exist separately in the Riemann invariants and v~ . 
(e) and (f) show the characteristics for f+ and v~ . The horizontal line at time 
t = 0.1 represents the remapping of the simulation. In both (e) and (f) the 
characteristics are converging to the shocks, while never intersecting. 



solutions of the homentropic Euler equations. It may be possible, however, 
to choose an averaging scheme such that the shock speeds between the two 
equations are the same or similar. 

Additionally the conditions under which a single traveling shock differs be- 
tween the two equations as is demonstrated in Examples 2 and 3. For a single 
shock to form with the CAHE equations either or v~ must be constant. 
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Even if it were possible to find an averaging scheme such that the shock speeds 
between the two equations are the same, these conditions would not change. 
For the homentropic Euler equations, the condition is tied in with the Rankine- 
Hugoniot conditions. 

The final difference noted between these equations is significant. As seen 
in Examples 2 and 3, for the homentropic Euler equations, there can be areas 
devoid of characteristics. These can be filled with expansion waves or by new 
characteristics originating from the shock. By averaging the characteristics in 
the CAHE equations, there will be no areas devoid of characteristics. Thus 
there will be no generation of new characteristics. This can cause significant 
differences in behavior as seen in Example 2, where v~ remains a constant for 
the CAHE equations, but not for the homentropic Euler equations. This again 
is a property that is unaffected by the averaging scheme chosen. 

8 Conclusions 

The CAHE equations were derived by taking the homentropic Euler equations 
and spatially averaging the characteristics. This led to a new set of equations 
that has many interesting properties. Existence and uniqueness proofs. Theorem 
|3.1| and [X2l were proven by establishing that the characteristics of the equations 
never intersect. In the process this established that any initial conditions in 
C''^(M) will have a solution in C^(M). Furthermore any discontinuities in the 
initial conditions will remain and be convected in the solution for all time. 

The speeds of shocks in the CAHE equations were found to be determined 
by the speed of the characteristics at the location of those shocks. Furthermore, 
different averaging schemes were shown to provide different shock speeds. 

The Riemann problem was then examined, where solutions were generally 
found to consist of two traveling discontinuities. However, often one of the 
discontinuities can prove to be unstable and by smoothing the initial conditions 
will develop into an expansion wave instead. 

Finally using some numerical examples and results from the previous sections 
the CAHE equations were compared and contrasted with the homentropic Euler 
equation from which they originated. Both can generate shocks and expansion 
waves from the Riemann problem. With the averaging scheme employed in the 
paper, the speed of the shocks differ. For the examples chosen, the behavior of 
the equations proved to be significantly different. Finally, the homentropic Euler 
equations showed the generation of new characteristics, while with the CAHE 
equations, new characteristics will never be generated. From this is seems clear 
that as a ^ the solutions to the CAHE equations will not converge to weak 
solutions of the homentropic Euler equations. 

The CAHE equations have proven to have convenient existence and unique- 
ness properties. When continuous initial conditions are chosen, the solution will 
remain continuous. However, the equations seem to display behavior with too 
significant of departure from the homentropic Euler equations to be of use in 
gas dynamic applications. It is possible that this could be rectified by taken a 
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one-sided average of the characteristics, but this concept was not pursued here. 
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